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This study evaluates the utility of various non-intrusive tech- 
niques for the measurement of the flow field on th^ windward side of 
the Space Shuttle or a similar re-entry vehicle. Included are 
linear (Rayleigh, Raman, Mie, Laser Doppler Velocimetry, Resonant 
Doppler Velocimetry) and nonlinear (Coherent Anti-Stokes Raman, 

Laser Induced Fluorescence) light scattering, electron beam fluores- 
cence, thermal emission and mass spectroscopy. Flow field proper- 
ties are taken from a nonequilibrium flow model by Shinn, Moss and 
Simmonds at NASA Langley. Conclusions are, when possible, based on 
quantitative scaling of known laboratory results to the conditions 
projected. Detailed discussion with researchers in the field con- 
tributed further to these conclusions and provided valuable insights 
regarding the experimental feasibility of each of the techniques. 
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NON- INTRUSIVE FLOW MEASUREMENTS 


ON A RE-ENTRY VEHICLE 

R. B. Miles, D. A, Santavlcca, and M. Zimmermar.n 

Abstract 

This study evaluates the utility of various non-intrusi ve techniques 
for the measurement of the flow field on the windward side of the Space 
Shuttle or a similar re-entry vehicle. Included are linear (Rayleigh, 
Raman, Mie, Laser Doppler Velocimetry, Resonant Doppler Velocimetry) and 
nonlinear (Coherent Anti-Stokes Raman, Laser Induced Fluorescence) light 
scattering, electron beam fluorescence, thermal emission and mass spectros- 
copy. Flow field properties are taken from a nonequilibrium flow model by 
Shinn, Moss and Simmonds at NASA Langley. Conclusions are, when possible, 
based on quantitative scaling of known laboratory results to the conditions 
projected. Detailed discussion with researchers in the field contributed 
further to these conclusions and provided valuable insights regarding the 
experimental feasibility of each of the techniques. 



NON- INTRUSIVE FLOW MEASUREMENTS 


ON A RE-ENTRY VEHICLE 


I . Introduction 

This project was undertaken to determine which if any of a variety of 
nonintrusive diagnostic techniques might be useful in determining flow pro- 
perties on the windward side of a re-entering Space Shuttle or a separate 
re-entry vehicle. Our procedure was to first estimate the detection capa- 
bilities of the techniques at conditions corresponding to those expected in 
the flow field surrounding the Space Shuttle as it passes between 80 and kO 
km in altitude. Where possible we used existing experimental results and 
extrapolated to the desired conditions. In some of these estimates we were 
guided by ongoing work in our own laboratories and in others by results pub- 
lished in the literature. Subsequently a presentation was made to the 
Princeton Advisory Group and the relative merits of each technique, plus a 
list of the important parameters to be measured were discussed. In order to 
further explore the promises and limitations of each of the techniques, a 
wide variety of researchers were visited and detailed discussions undertaken. 
These inclu >ed trips to NASA Langley, United Technologies, NASA Ames, Uni- 
versity of Cal i fornia-Berkeley, Sandia-Livermore Laboratories, and Stanford 
University. This provided us with an opportunity to learn the state-of-the- 
art in each of the areas and to obtain specific comments and perceptions 
from these researchers with regard to the potential of each of the techniques. 
Where relevant many of these comments are included in the appendices in 
which the techniques are discussed in detail. 
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2. Summary of the Results 

The flow field conditions which were used in this study were largely 
taken from work conducted at NASA Langley by J. L. Shinn, J. N. Moss and 
A. L. Simmonds.' Figures 1 and 2 show plots of velocity, density and 
temperature for a re-entering Shuttle at altitudes of 7^ .95 km and 52 .A3 
km, respectively. Each is shown as a function of the distance from the 
wind-ward surface at two points along the centerline. Here, the calcula- 
tion is done for an equivalent hyperboloid whose parameters are determined 
by the pitch angle of the re-entering vehicle. The arrows indicate the 
conditions at the wall assuming a finite catalytic activity. The free 

2 

stream conditions and equivalent body parameters are given in Table 1. 

Each of the curves in Figures 1 and 2 terminates at the shock boundary and 
thereafter the free stream conditions apply. 

These calculated results depend strongly on the effects of nonequilib- 
rium heating and surface catalytic activity. For example, Figures 3 and 4 
show calcu’ated values of species mass fractions at 75 km altitude assuming 
in Figure 3 a non-catal y ti c wall and in Figure 4 a fully or equilibrium 
catalytic wall. Depending on the range in catalytic activity, the wall 
temperature can vary by over 30%. In Figures 3 and 4 note that there is a 
large change in the species concentration near the wall with a large NO 
concentration indicating strong catalytic effects. At low altitudes the 
temperature behind the shock is too low to produce significant radical 
concentrations (Figure 5).^ These species and temperature distributions 
behind the shock depend strongly on estimates of reaction rates. Conse- 
quently, the validity of the model remains largely unknown until temperature 
profiles or species profiles can be measured. Density and velocity profiles 
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Figure 2. Velocity, temperature and density outward from the windward 

centerline of the re-entering Space Shuttle at an altitude of 
52.1*3 km. Conditions are shown for two locations (s*) : 1.1 

meters from the nose (s/R^- .9) and 30.8 meters from the nose 

(s/R n « 25) [Ref. 1]. 
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Figure 4. Species mass fraction for nonequilibrium flow with an equilib- 
rium catalytic wal I [Ref. 1]. (1*2. 2° Hyperboloid; R*. ™ 1 . 3A 1 m; 

j/R n - 10.0; Alt- 7**98; M^-27.5; T w -1500°K). 
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Figure 5. Species mass fraction for nonequi Hbrl urn flow with a finite 
catalyti c wal 1 [Ref. 3l» Alt* 156.5 ft (^7*67 km); M*9«15» 
a- 3**. 8°; s/R N * 2 5* 
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Table 1. 

Freest ream 

conditions and 

equivalent body 

paramete rs 

L 

Case 

Altitude 

km 

• u > 

1 00 ’ 

km/sec 

Poo* 

kg/m 3 

H eo 

a 

Equivalent Body* 
e r n , m 

1 

74.98 

7.62 

3.80x10 

27.5 

41.4 

42.20 

1 .341 

2 

68.88 

6.61 

9.28xlO -5 

22.6 

40.2 

40.85 

1.286 

3 

64.92 

5.94 

1 .55x10 _i * 

19-6 

34.5 

35.20 

1 .042 

1* 

60.66 

5.12 

2.6lxl0 _i * 

16.4 

31.8 

32.40 

.938 

STS-l Conditions 

5 

52.43 

3.62 

7.43xlO -/ * 

11.1 

39.0 

39.7 

1.233 

6 

48.16 

2.98 

1 .33xlO -3 

9-15 

34.8 

35.4 

1.052 

7 

44.04 

2.30 

2.4Cxl0‘ 3 

7.14 

30.9 

31 .6 

.908 


a angle of attack 

* Hyperboloid with body half angle 8 and nose radius at zero degree 
angle of attack. 
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as well as the location of the shock are also of significant interest. Near 
the nose the shock stands less than 20 cm from the surface and substantial 
variations in velocity, temperature, density and species concentrations oc- 
cur over a few millimeters. These gradients are largest through the boun- 
dary layer where phenomena most directly related to the surface will occur. 

The capabilities of the various nonintrusive detection methods are 
shown in Table 2. This table breaks the techniques into three categories. 

The first includes passive techniques which do not require illumination of 
the flow. Thermal emission naturally occurs from hot gases aid particles 
and may be collected through a window and analyzed by an optical spectrometer. 
Hass spectroscopy involves sampling at the surface and subsequent analysis 
of gas constituents. The second category includes techniques which require 
actively probing the flow field with light or an electron beam. The last 
category includes techniques which require both active probes and a foreign 
species seeded into the flow field. Note that, in theory, numerous ways are 
available for measuring both temperature and species concentrations. 

The selection of the flow properties which would be the most valuable 
to measure was a topic of substantial discussion. Clearly one would like 
to measure everything. Almost any measurement would give some check on the 
accuracy of the modelling. It was generally agreed however that species 
concentration or temperature would be especially useful as a check of the 
model. Certainly a profile of one of these two variables would produce data 
of substantial interest. This is not to neglect the importance of the den- 
sity profile, shock location, velocity profile or the measurement of free 
stream conditions. One of the most common concerns raised was the density 
of particles and the amount of thermal emission encountered in the flow. 



Table 2 

Capabilities of Nonintrusive Detection Methods 


TECHNIQUE 

MEASUREMENT 




1 . Mass Spectroscopy 


Spec 1 es 
Concentration 



Thermal Emission 

Tempe rature 

Species 
1 denti ty 



II. Rayleigh Scattering 

Temperature 


Density Velocity 


Raman Scattering 

Tempe rature 

Species 

Concentration 



Coherent Anti -Stokes 
Raman Scatter! ng - CARS 

Tempe rature 

Species 

Concentration 

Den si ty 


Mie Scattering 




Particle 

Size/Density 

Laser Induced Fluorescence 

Temperature 

Species 

Concentration 



Electron Beam Fluorescence 

Temperature 

Species 

Concentration 




Laser Doppler Veloclmetry 


Veloci ty 

Particle 

Size/Density 

Resonant Doppler 
Veloclmetry 

Tempe rature 

Density Velocity 
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These two questions clearly affect the viability of several of the tech- 
niques and should be answered at an early stage. As will be evident In the 
following discussion, species measurements will be the most difficult to 
obtain with existing techniques. On the other hand, density profiles, 
temperature profiles, velocity profiles, particulate profiles, and the 
shock location may well be achievable with current technology. 

A detailed discussion of each of the techniques and the methodology used 
for evaluation is presented at the end of this report. Table 3 is a summary 
of the conclusions. This table Is broken into five categories. The first 
includes techniques which have undergone substantial laboratory testing and 
appear to have high promise. These are Rayleigh scattering, especially for 
density measurements and potentially for temperature and velocity measure- 
ments, and Mie scattering for particle concentrations, velocity and size. 

The second category includes techniques which also appear promising but have 
either not been adequately demonstrated in the laboratory or their success 
depends strongly on the magnitude of the background noise levels. These are 
laser induced fluorescence measurements of N, 0 and NO concentrations and 
temperature (from NO), and mass spectroscopy. Laser induced fluorescence 
may produce profiles of the N, 0 or NO concentrations which are sensitive 
indicators of the chemistry at high altitudes. Hass spectroscopy may be 
capable of determining species concentrations within the boundary layer 
throughout re-entry. The third category includes techniques which will 
require substantial additional work and are of questionable utility. These 
are quantitative thermal emission, electron beam fluorescence, and spontan- 
eous Raman scattering measurements of temperature and concentration. The 
fourth category includes techniques which appear to have very little promise. 
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Table 3. 

Summary of Conclusions 


CATEGORY TECHNIQUE 


1 

A. 

Rayleigh Scattering 


B. 

Hie Scattering 

2 

C. 

Laser Induced Fluorescence 


D. 

Hass Spectroscopy 

3 

E. 

Thermal Emission 


F. 

Electron Beam Fluorescence 


G. 

Raman Scattering 

A 

H. 

Coherent Anti -Stokes Raman 
Scattering - CARS 

5 

1. 

Laser Doppler Ve loci met ry 


J. 

Resonant Doppler Velocimetry 


PROJECTION 

Very Promising 
Very Promising 

Promising 

Promising 

Needs More Research 
Needs More Research 
Needs More Research 

Little Promise 

Require Seeding 
Require Seeding 



- 13 - 


The only technique in this category is coherent anti-Stokes Raman scattering 
where the low densities, high temperature and crossed-beam phase-matching 
result in low signal strengths. A fifth category is reserved for those 
techniques which require seeding. It is currently impossible to determine 
the promise of these techniques, however it may be possible to use ablated 
material or naturally occurring particles in such a manner that the tech- 
niques may be viable. These are laser Doppler velocimetry and resonant 
Doppler velocimetry. Each of the techniques is further discussed in its 
associated appendix. 

In our estimation the most promising technique is Rayleigh scattering. 
This has been demonstrated in the laboratory to be capable of measuring 
density and temperature. The scattering is observed at the same frequency 
as the laser and therefore the signal is subject to substantial noise inter- 
ference from scattering from particles which may be in the flow stream. 
Similarly, observations very close to the window may be difficult due to 
background light scattered from the window itself. Various species have 
different Rayleigh scattering cross sections. These do not vary substan- 
tially from that of and, since Nj is always the major species, its cross 
section may be used for density and velocity measurements. For temperature 
measurements, however, the different molecular weights of the different 
species will cause different broadening of the scattered light, somewhat 
complicating the measurement. The Rayleigh technique involves a single 
laser probe and could use an array of detectors so that a single laser pulse 
would yield an ertire profile. It has been demonstrated in the laboratory 
that such a system may be integrated with a Mie scattering system for the 
detection of particulates. For example, when light is strongly scattered 
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from a particle, simple logic in the detection system may be used to exclude 
that signal from the Rayleigh scattering measurement and separately register 
it in a particle counter. Thus a single apparatus may conceivably yield 
pressure, temperature, velocity and particulate profiles. 

We also find laser induced fluorescence from 0, N and KC to be a prom- 
ising technique for both concentration and temperature measurements. This 
has not yet been fully demonstrated in the laboratory under conditions which 
can easily be related to the re-entry environment. Therefore we must in- 
clude this technique in the second category. Our enthusiasm is based on the 
importance of these species as an indicator of the nonequilibrium chemical 
dynamics, particularly at high altitude. Although there are difficulties in 
using this technique to measure species concentrations, the wide range of 
species densities behind the shock should nonetheless be observable. The 
work of McKenzie et al. at NASA Ames has been largely directed towards the 
use of this technique for temperature measurements. Using the technique as 
he has done, it would appear that temperature profiles may be obtainable 
using two sequential laser pulses of different wavelengths. With this tech- 
nique we also note that scattering may be observed by multiple detectors as 
was true of Rayleigh scattering so that a profile may simultaneously be 
taken. Since a pulsed laser is used and fluorescence is observed over a 
rather narrow range of the spectrum, thermal emission can largely be excluded. 

Our evaluation of mass spectroscopy was limited since no operational 
system for extracting the gas sample through the vehicle surface has yet 
been developed. There are two serious limitations associated with this 
technique. The first is that it samples from an area in such close prox- 
imity to the sampling port that it may not be Indicative of the conditions 
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at the surface at all. This will particularly be true if the sampling port 
is large in the streamwise direction. The second problem Is that the sample 
gas in passing through the sampling probe will most certainly undergo recom- 
bination and as such will not be truly indicative of the species concentra- 
tions in the boundary layer. Discussions with Dr. William McLean at Sandia 
and Drs. Daniel Seery and Med Colket of United Technologies indicate that a 
scavenging probe may be useful to alleviate this problem. For example, one 
could introduct hydrogen at the sampling port to interact with the oxygen 
atoms or Ojg to interact with nitrogen atoms to produce a stable species that 
would in effect tag the original species. A similar idea is being studied 

by George Wood at NASA Langley using isotopic labelling on the wall material. 

4 

This technique was originally suggested by Fristrom and is currently being 
used for radical species diagnostics in combustion. With proper care in the 
design of the intake port and a careful selection of pumping rates, this may 
be an effective way of doing species identification in close proximity to 
the surface. 

Although thermal emission and electron beam fluorescence may have the 
capability of yielding temperature and species concentrations, they are both 
better used for qualitative than quantitative study of the flow field. For 
example thermal emissions might enable one to determine the total column 
density which could then be used with laser induced fluorescence to correct 
for I'-apping effects as has been suggested by Professor John Daily of the 
University of California. The thermally emitted light must be collected 
over an integrated path and therefore point measurements will be difficult 
if not impossible. Pressure and temperature effects will further confuse 
the analysis. Similarly electron beam fluorescence has been developed as a 
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flow visualization tool and is not well adapted to quantitative measurements 
under these conditions. As Dr. William Hunter of NASA Langley pointed out, 
both quenching and temperature effects are not well understood. The elec- 
tron beam does have the advantage however that it will operate at very high 
altitudes and might be useful as a qualitative density probe for the low 
densities encountered in the early stages of re-entry. 

The quantum noise limited sensitivity estimates for spontaneous Raman 
predict that N, 0, NO, 0^ and N^ concentration measurements and N^ vibra- 
tional temperature measurements should be possible over the entire range of 
altitudes. This is encouraging but perhaps misleading since the smallest 
amount of background noise from thermal emission or laser induced particu- 
late incandescence will render this technique useless. Until more informa- 
tion is available about the thermal emission and particle density conditions 
however it is reasonable to keep Raman scattering under consideration. 

Coherent anti-Stokes Raman scattering was found to be incapable of 
providing useful information under these conditions. The combination of low 
density, high temperature and short interaction length (because of the 
required crossed beam phase matching) lead to a signal strength below the 
detection limit. It is also restricted to single point measurements and, 
because of the required wide angle crossed beam phase matching, it would be 
very difficult to measure profiles unless a tremendous amount of experimen- 
tal complexi ty were involved. 

If the flow field could be seeded either from the surface of the vehi- 
cle or by a cloud of seed material, then several other techniques become 
available. In particular, laser doppler velocimetry (LDV) may be an attrac- 
tive way of measuring velocity profiles. Due to the very high velocities, 



however, conventional LDV configurations are not appropriate. Several al- 
ternatives may be considered including a simple monitoring of the transit 
time scattering intensity profile as a particle passes tnrough a single 
light beam as has been discussed by Dr. Don Holve of Sandia-Livermore, or the 
use of interferometers to observe the large spectral shifts off the light 
scattered by rapidly moving particles. Another possibility is the crossing 
of two collimated laser beams at a very slight angle to produce widely spaced 
interference fringes to lower the modulation frequency of the fast particles 
passing through the laser beams. In this case the depth of field becomes 
large and a multiple detector system may be able to observe the velocity 
prof i le. 

Similarly if sodium atoms could be seeded into the flow, then the reso- 
nant Doppler velocimeter could be used to measure temperature, density and 
velocity. This technique requires that the laser frequency be scanned and 
thus integration times of e few seconds would be expected. Here again an 
array of detectors can simultaneously observe the entire flow field profile. 
Using this technique the velocity profiles should be easily discernible, 
however quenching effects and the wide range of temperatures may lead to 
difficulty with the temperature and density measurements. 

3. Recommendations 

Throughout this study many questions were raised concerning the density 
of particulates and the amount of thermal emission on the windward side of 
the re-entry vehicle. It should be clear that the density of particulates 
is au extremely important quantity since it provides both a noise source and 
a potential seed source for laser Doppler velocimetry. In the case of Raman 
scattering even relatively low levels of thermal emission or low densities 
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of particles will prohibit Its use. Clearly If the number of particulates 
Is exceptionally high, it may also obscure the Rayletgh scattering or the 
laser induced fluorescence signals and eliminate laser Doppler velocimetry 
as an effective diagnostic tool. On the other hand, if the number of parti- 
culates is low, Rayleigh scattering may be viable and might be used for 
temperature and velocity as well as density measurements, but laser Doppler 
velocimetry may not be possible. Similarly, thermal emission, if excessive, 
may obscure the laser induced fluorescence signal as well as Rayleigh scat- 
tering. For these reasons we recommend that a preliminary experiment be 
planned to measure the particle densi* and therr.-,l emission. This could be 
done using a simple laser Learn and through a window collecting the Mie scat- 
tering at the laser frequency and the broad band thermal emission. As pre- 
viously mentioned, this could be integrated with a Rayleigh scattering exper- 
iment if proper signal processing was done. With this information, prelim- 
inary estimates of the particle density and thermal emission intensity could 
be made and perhaps a density profile found. 

Our results indicate that Rayleigh and Mie scattering and laser induced 
fluorescence are promising techniques. We would like to recommend that fur- 
ther work be undertaken on the development of these techniques with this 
particular application in mind. With regard to mass spectroscopy, our con- 
clusion is that the technique has promise but is as yet unproven. It would 
appear that a measurement of species concentrations would require the use of 
a reacting tag, either in the gas phase, or on the surface of the probe. 

The major limitation of this technique Is its inability to probe beyond the 
surface and the possibility that the catalytic nature of the entrance port 
may perturb the species concentrations from what would normally be encoun- 
tered along the vehicle surface. 
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APPENDIX A. RAYLEIGH SCATTERING 

Rayleigh scattering is an elastic scattering process for which the scat- 
tering intensity is proportional to the number density of scattering moler 
cules and therefore can be used to measure density. It is several orders of 
magnitude stronger than Raman scattering but lacks the species selectivity 
of Raman scattering and is easily obstructed by Hie scattering (in particle 
laden environments) and reflected laser light. It has been used for temper- 
ature measurements in constant pressure, non-reactive flows where the densi- 
ty and temperature are related through the equation of state [A.l]. The 
Doppler shift and broadening of the Rayleigh line can also be used to mea- 
sure the velocity and temperature, respectively [A. 2]. The Rayleigh scat- 
tering cross-section is species dependent, therefore a change in Rayleigh 
signal strength or line shape can be a result of composition changes as well 
as density or temperature changes. 

The Rayleigh scattered signal strength, in terms of power, is given by 
the following expression: 

P RAYLEIGH “ P o n C n D | N t Q i 

where P Q is the power of the Incident laser beam, L Is the length of the in- 
cident laser beam from which scattered light is collected, fi c Is the solid 
angle of the collection optics, is the efficiency of the collection optics, 
Hq is the detector quantum efficiency, Nj is the number density of species 
i and Qj is the Rayleigh cross section for species i. 

The sensitivity and accuracy of Rayleigh scattering for total density 
measurements was evaluated in terms of the minimum detectable density for a 
given signal to noise ratio. The noise was assumed to be quantum noise 




A- 2 


which Is equal to the square root of the number of detected Rayleigh photons. 
This leads to a signal to noise ratio: 


(S/N) 


RAYLEIGH 


^RAYLEIGH At n p /hv * 


1/2 


(A. 2) 


where At Is the pulse duration, n Is the number of pulses during the mea- 

P 

surement interval, h is Planck's constant and v is the radiation frequency. 

The experimental parameters used in equations (A. I) and (A. 2) to deter- 
mine the minimum detectable total density are those given In Table A.) Two 
lasers were considered, a krypton fluoride (KrF) exclmer laser and a copper 
vapor laser. The KrF because of its high average power and its UV output, 
which results in larger cross-sections and greater detector quantum effi*' 
ciencies. The copper vapor because of its narrower output linewidth and 
higher efficiency. Based on the rate of descent of the Shuttle, an integra- 
tion time of 10 seconds was chosen. The Rayleigh cross section for was 

_2 p 2 

used in this estimate, l.e. 8x10 J cm /str-molecule at 5l*t.5 nm. For a 
signal to noise ratio of 10, the. resultant minimum detectable density is 
3x10 /cc with the KrF laser and 6x10 /cc with the copper vapor laser. 
Both of these values are below the minimum density at the highest altitude. 
This result agrees well with a detection limit estimate made by Robert 
Dibble of Sandia-Ll vermore Laboratory who extrapolated from detection 
limits he routinely observes In the laboratory. He also estimated the 
error due to variations in composition to be less than ± 1 5% . Therefore 
assuming we are in the quantum noise limit the use of Rayleigh scattering 
for total density measurements over the entire range of altitudes appears 
very promising. The problem of Hie scattering from particulates Is an 
important concern since it will be orders of magnitude greater than the 



Table A.1 


Values Used for Rayleigh Scattering Signal to Noise Estimates 




KrF Laser 

Copper Vapor Laser 

Peak Power 

P 

o 

40 MW 

0.33 HU 

Pulse Length 

At 

25 nsec 

30 nsec 

Repetition Rate 

n 

P 

l0 2 /sec 

5 x 10* /sec 

Wavelength 

A 

249 nm 

5)0 nm 

Detector Quantum Efficiency 

*0 

0.35 

0.1 

Collector Optics Efficiency 

^C 

0.1 

0.1 

Optics Collection Angle 

ft c 

0.01 str 

0.01 str 

Scattering Length 

L 

5 mm 

5 mm 



Rayleigh scattering; however, if the particle density is not too high it is 
possible to identify and reject the larger signals due to particles [A. 3]. 

To use Rayleigh scattering for temperature and velocity measurements, 
it is necessary to use a narrow line width laser (which rules out the use 
of the KrF laser) and a scanning Fabry-Perot spectrometer with adequate 
resolution, in order to resolve the Doppler shift and broadened line width 
which will both be of the order of 0.1 cm The number of detected photons 
will in turn be reduced, which will increase the minimum detection limit to 
approximately lO*^/cc. Since the minimum nitrogen density at 72 km altitude 
is also 10^/cc, the Rayleigh temperature and velocity measurements will be 
possible at all altitudes, throughout the flow field. For the temperature 
measurement it will be necessary to correct for composition changes perhaps 
using the calculated Rayleigh lineshape [A. 2], although the densities are so 
low that this latter effect is likely to be negligible. 
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APPENDIX B. HIE SCATTERING 

The presence of particles which have ablated from the space shuttle 
tiles are of Interest for several reasons. One is that they represent a 
potential noise source for all of the optica! diagnostics. Mle scattering, 
together with laser induced particulate incandescence and thermal emission 
could interfere with other optical measurements depending on the particle 
density and size distribution. It is essential to obtain information about 
the severity of these potential problems before implementing any of the op- 
tical techniques. This would require installing a window in the windward 
side of the re-entry vehicle. This window could be used to measure the 
intensity and spatial characteristics of the thermal emission from both the 
gas and the particles, as well as to look at the Hie scattering from a laser 
to get an estimate of the particle density, and, possibly, the size and spatial 
distribution. 

The particles can also be used as a diagnostic seed for measuring velo- 
city and temperature. Shortly after the particles are ablated from the tile 
surface they will be accelerated to the gas velocity. The time required to 
reach the gas velocity will depend on the particle mass and the gas velocity. 
The depth of penetration into the boundary layer will depend on the perpen- 
dicular momentum with which the particle was released from the surface. 

Within a certain distance from the surface there will be some particles 
which have not yet equilibrated with the gas velocity. It is necessary to 
determine which particle velocities correspond to the gas velocity. One 
approach is to simply take the velocity of the fastest particles to be that 
of the gas. Another method is to measure both the size and the velocity 
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simultaneously and to assume that the small particles are most likely to be 
travelling at the gas velocity. Such Information could be obtained from a 
dual beam laser Doppler velocimetry system where the particle size is ob- 
tained from the visibility parameter of the Doppler burst. Unfortunately 
this method requires a priori knowledge of the particle size range since 
the visibility parameter is single valued over only a limited range of 
particle sizes. No other simultaneous velocity and particle sizing tech- 
nique has been demonstrated to a reasonable degree of satisfaction. Dr. 
Donald Holve of Sandia-Li vermore suggested using single beam methods, l.e. 
transit time velocimetry for velocity and pulse height detection for parti- 
cle size, whereby velocity and size distributions could be obtained but not 
velocity-size correlations [B.1,B.2]. Techniques where one looks at the 
angular distribution of the Hie scattered light would be difficult to inter 
pret because of the unknown and irregular particle shapes. Additional de- 
tails on velocity neasurement techniques can be found in the sections on 
laser Doppler velocimetry and resonant Doppler velocimetry. Another possi- 
bility for measuring the particle velocity which was suggested by Dr. 

Holve is to use double pulse or streak photography where, for example,a 
particle travelling at I km/sec would travel I mm during a 1 microsecond 
light pulse. Spatial resolution would be determined by the depth of field 
of the camera optics. 

The particles could also be used to give an indication of the gas tem- 
perature by looking at the spectral distribution of the thermal emission 
from the particles. One would have to make some rssumptions about the par- 
ticle emissivity to determine its temperature and then the question of ther 
mat equilibrium between the particle and the gas would still have to be 
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considered. As discussed in the section on Rayleigh scattering, a combined 
Hie scattering and Rayleigh scattering experiment looks promising as a 
method for determining particle densities. 
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APPENDIX C. LASER INDUCED FLUORESCENCE 


Laser Induced fluorescence has the potential of measuring rotational 
and vibrational temperature and species concentration and of performing 
density sensitive flow visualization. Generally one '-r more laser beams 
are directed into the flow to excite a particular species to a higher ener- 
gy level. A detector then collects the fluorescence as the gas relaxes 
back to its initial state after a characteristic lifetime. The detector may 
be configured to monitor any position along the laser beam to give good spa- 
tial resolution. The fluorescence intensity is proportional to the popula- 
tion of absorbing molecules in the radiating upper energy level during the 
excitation process. This population is related to the total population 
through the rate equations which include the processes of spontaneous emis- 
sion, stimulated emission, absorption, collision quenching, rotational re- 
distribution, and vibrational redistribution. Unfortunately, the parameters 
describing these competing processes are not always known; particularly, 
the identity, concentration, and quenching rates of all possible collision 
partners are important and largely unknown quantities. 

Temperature and species concentration measurements using laser induced 
fluorescence have been demonstrated at pressures equal to or less than one 
atmosphere and at temperatures up to those found in typical combustion 
environments (about 2000°K) . The most frequently treated species has been 
OH, because of its important role in combustion kinetics; however, work has 
also been performed on numerous other systems including NO, N and 0 [C.l - 
C.*»] which are of interest in this study. 

Several schemes have been proposed to account for the effect of 
quenching on the experimental results. One method, as first proposed by 
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Piepmeier [C. 51, requi res saturating the upper energy level and thereby 
eliminating the dependence on laser Intensity as well as the need to cor- 
rect for quenching. Since complete saturation is difficult to achieve, 
Baronovski and McDonald [C.6] proposed and demonstrated a partial saturation 
technique whereby one can extrapolate to the complete saturation limit. 
Another method, proposed and demonstrated by Sterowski and Cottereau [C . 7] 
does not require saturation but requires that the laser pulse duration be 
shorter than the time for steady state conditions to be reached. The quench 
ing rate then can be related to the fluorescence decay curve. The required 
laser pulse duration and the fluorescence decay rate will depend on the gas 
pressure and temperature. Other methods include the use of empirical quench 
ing corrections as in the work of Muller et al . [C.8], and the use of de- 
tailed multi-level models as in the work of Chan and Daily [C . 9 ] and that 
of Kotlar, Gelb and Crosley [C . 10]- 

The most promising approach given the low density conditions is to 
determine the quenching rate from the fluorescence decay rate as has been 
done by Stepowski and Cottereau [C.7] for OH in a low pressure flame and by 
Bischel, Perry and Crosley [C.ll] for both 0 and N atoms in a low pressure 
flow discharge. Based on the results of the latter work and the similarity 
between the excited electronic state lifetimes for OH and NO, a laser pulse 
width of less than 5 nsec is required. For 0 and N atom fluorescence, 
two-photon absorption is required because the first resonance transition 
lies far into the VUV. As Bischel et al [C.ll] have demonstrated, one 
would detect the resultant near IR fluorescence. For NO the r esonance 
transition can be accessed by either single or two-photon absorption. The 
disadvantage of the two photon process is that it is non-linear in laser 
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power. The advantage is that it is less difficult to generate tunable 
stable, high power visible radiation than the UV radiation required for the 
single photon process. There was no concensus among the researchers inter- 
viewed on this point, in either case the best laser system would be a 
Nd:YAG pumped dye laser using frequency doubting and/or Raman shifting to 
generate the required UV wavelengths. 

After quenching has been accounted for, the fluorescence intensity 
gives the population of the lower energy level before laser excitation 
which can then be related to the total population if the temperature is 
known. The fluorescence technique can also be used to find the tempera- 
ture by measuring the population of two or more energy levels since the 
ratio of the population of any two levels is equal to the ratio of the 
corresponding Boltzmann factors. 

Of particular interest to this study is the work of McKenzie and Gross 
[C.2] who have demonstrated single pulse temperature measurements of NO 
seeded into a flow of nitrogen. A selected rotational level is two photon 
excited and its fluorescence monitored. A very short time thereafter a 
second laser pulse two photon excites the NO from a different rotational 
level to the same upper state. The ensuing fluorescence intensity due to 
each pulse can be related to the rotational populations in the two ground 
states, the ratio of which can be related to the rotational temperature. 

Since only the ratio of energies absorbed needs to be known, the redistribu- 
tion of energy among rotational and vibrational levels of the excited state 
and their non-radiative quenching has a negligible effect on the measurement. 
This is true provided the quenching rate is insensitive to the rotational 
quantum number and broad-band detection is used. McKenzie considered the 



transition A^E + , v'-O X^TI ^ * v"»0 in NO from the J" “ rotational 

level since it has the largest absorption coefficient at room temperature 

and is separated far enough from other lines to allow selective excitation. 

He obtained a signal to noise ratio (SNR) in the quantum noise limit of 200 

for a single 1 mil li joule pulse energy of 5 nsec duration. The NO concen- 
ts -3 

tration which he used was 2x10 cm at room temperature and the size of 

the probe volume was 0.1 mmxl.O mm. The maximum density of NO in the flow 
surrounding the re-entry vehicle is an order of magnitude higher at an alti- 
tude of 75 km. At this higher temperature the NO molecules distribute over 
more states thus lowering the density in any given vibrational rotational 
level. The signal is further reduced because of the increased quenching due 
to the higher temperature. 

IS -3 

The minimum NO density downstream of the shock is 10 cm and the NO 
may have a translational temperature as high as 16,000°K. Scaling calcula- 
tions were made for both 16,000°K and 3000°K and yielded a SNR of 9 and SO 
respectively for a single 1 mJ pulse. The lower temperature was chosen due 
to the expectation that internal modes and chemical reactions will rapidly 
cool the translational energy of the gas, substantially lowering the rota- 
tional temperature. Increasing the laser's intensity to improve the SNR is 
limited by Stark broadening. The broadening reduces the laser energy coup- 
ling into the transition and may cause several lines to coalesce so that a 
particular rotational level may no longer be preferentially pumped. 

McKenzie has not yet quantified this limitation. 

To determine the density detectability limit for a quantum noise limited 
signal to noise ratio of 10, a laser operating at 10 pulses per second with 

1 mJ per pulse was considered. The integration time was assumed to be 10 

IT -T 

seconds. The NO detectability is then found to be 1.5x10 J cm J and 7.5 x 



1 1 - ? ? 

10 cm J for a temperature of 16x10 K and 3000°K respectively. The NO 
number density in most of the region between the shock and the vehicle sur~ 
face at an altitude of 75 km is predicted to be between 10^ cm ^ to 10^ 
cm . A comparison of Figures 4 and 5 in the main text of this report 
shows that the NO concentration drops quickly as the altitude decreases and 
by 47 km it is negligible except in the boundary layer. Consequently the 
laser induced fluorescence technique is only useful at the higher altitudes. 

In summary measurements of NO, N and 0 concentrations and NO rotational 
temperature appear feasible at 75 km where appreciable concentrations are 
expected. Such measurements have been demonstrated under similar but not 
identical conditions. Thermal emission from the high temperature gas behind 
the shock and hot particles as well as Hie scattering from particles should 
be investigated as a potential problem. It would therefore be advisable to 
monitor the environmental conditions before a laser induced fluorescence 


measurement is carried out. 
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APPENDIX D. MASS SPECTROSCOPY 


Hass spectroscopy can be used to separately measure the concentrations 
of 0 2> N NO, 0 and N over the entire range of ‘ensitles that are expected 
during re-entry. In addition, using a mass spectrometer with a fixed col- 
lector plate detector, as has been proposed by George Wood of NASA Langley, 
would enable simultaneous measurements of all these species. A similar 
spectrometer has been built and successfully tested by Lehotsky [D.l]. Dr. 
Martin Zabielski of United Technologies pointed out that the sample flow 
rate between the altitudes of 80 and <>0 km would change by approximately a 
factor of 10 for a fixed orifice area and given pumping system. He suggested 
that the pumping system be designed for kO km and then simply overpumped at 
higher altitudes. This would permit continuous measurements to be made at 
all altitudes. 

If effusive collection is used, (the collection orifice diameter less 
than the mean free path) then the sample would come from within a few mean 
free paths of the surface. This may cause difficulty since the gas composi- 
tion within a few mean free paths of the collection orifice may be affected 
by the presence of the orifice which Is likely to be made of a different 
material and at a different temperature than the surrounding tiles. Dr. 

David Seery of United Technologies suggested that it may be possible to use 
a larger orifice and collect the sample by modest pumping which, because of 
the low densities, would cause minimal distrubance to the flow field. By 
pumping on the orifice one would extend the sample volume further out into 
the flow field, beyond the region which is affected by the local surface 


reactions. 
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A major problem In using mass spectroscopy is the extraction of a sam- 
ple which reflects the true composition of these gases. Radical species, 
particularly 0 and N atoms, quickly react or recombine. Because of the 
thckness of the tiles protecting the surface of the re-entry vehicle and 
the necessity to thermally insulate the mass spectrometer from the exterior 
temperatures, the sample must travel several inches before being analyzed. 

In order to minimize the effect on the local flow field, the passageway 
through the tiles should be small, e.g. a differential pumping, skimmer 
collection system could not be used. With the sample having to travel 
through a long, narrow passageway, both wail and gas phase collisions, and 
in turn 0 and N atom recombination reactions, will be inevitable. Wood has 
proposed selecting the wall material such that the 0 and N atoms selectively 
react upon collision with the wail producing a stable and identifiable spe- 
cies which would in effect serve to tag the original atomic species. Drs. 
Daniel Seery and Med Colket of United Technologies and Dr. Robert McLean 
of Sandia-Li vermore , all suggested using the scavenging piobe concept pro- 
posed by Fristro;n and recently demonstrated by Fristrom and McLean [D.2]. 
Here, the sample is immediately diluted at the orifice with a gas stream 
containing a species which selectively reacts with either the 0 or N atoms 
to produce a stable and identifiable species. This species then serves to 
tag the original atomic species. If this stable product were already pre- 
sent in the sample then it might be possible to distinguish between the 
sample and the tag molecule by using different isotopes. 

Another possibility suggested by Seery and Colket Is to extend a tra- 
versing scavenging probe out into the boundary layer in order to extract 
samples at different distances from the wall. They felt that because of the 
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very high velocities that all appreciable disturbances to the flow field 
would be convected downstream and would not affect the measurement. They 
did point out, however, that when the flow is supersonic, as It Is outside 
the subsonic boundary layer, there will be a shock in front of the probe 
which may substantially alter the sample's composition. 

The general conclusion of this review is that using a scavenging probe 
to extract the sample with subsequent mass spectroscopic analysis Is a pro- 
mising technique for determining the concentrations of N^, 0^, 0 and N near 
the wall and even possibly somewhat away from the wall into the boundary 
layer by using a traversing probe or by modest pumping rates. 
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APPENDIX E. THERMAL EMISSION 


Thermal emission from atoms and molecules may be used to measure tem- 
perature, identify species and perform density sensitive flow visualization. 
Thermal emission from hot particles may also be used to obtain the particle's 
temperature and velocity. It may alsc constitute a noise source for other 
optical nonintrusive techniques. 

The fact that the collected signal is integrated over a path length 
prevents this technique from yielding good spatial resolution. Since most 
flows are inhomogeneous, an integration must be performed on the flow param- 
eters as well. The problem is compounded if optically thin and optically 
thick regions are present along the path since the method for analyses of 
each of these regions is different. A major source of uncertainty in the 
optically thin limit is due to incomplete knowledge of quenching rates and 
Franck-Condon factors. 

In a nonequilibrium gas, the concept of temperature is not defined. 
Usually the translational modes quickly equilibrate and thus most accurately 
reflect the temperature of the local environment. The translational temper- 
ature of an optically thin gas can be measured by resolving the spectral line 
profile of the emission from a particular transition. The profile has to 
be taken by using an interferometer to achieve high enough resolution. For 
a purely thermally broadened line, the frequency full width at half maximum, 
Av, is related to the gas temperature, T, by [E.l], 


T 


, Av A 2K 2 H 
1 2 ’ 2k ln2 


(E.l) 


where and M are the emission wavelength and mass of the particular spe- 
cies and k is Boltzmann's constant. When using this relation to measure T 
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it is essential to account for the instrumental width and collision broad' 
ening. 

The population of the various energy levels may reach separate thermo- 
dynamic equilibria which differ from the translational temperature and 
persist for a relatively long time. The temperature parameter in the Boltz- 
mann distribution of energy level populations may be measured in a variety 
of ways. The emission intensity, from an optically thin source using 

the Boltzmann population factor is [E.2] 


*21 ' K 9 2 '*P<-E 2 ' kT) 


(E.2) 


where K is a constant, g2 is the degeneracy of the excited level, is the 
Einstein spontaneous emission coefficient of the transition, and Ej is the 
energy of the excited state. Apart from uncertainties regarding the opti- 
cal depth, the major errors involved in using this method are uncertainties 
in the Einstein A coefficients and quenching rates. Equation (E.2) may be 
applied to electronic, vibrational, or rotational population levels depend- 
ing on the excitation energy. 

Both the translation and the population temperature measurement tech- 
niques depend on the gas being optically thin. As the gas turns optically 
thicker, the emission spectrum approaches the black body curve as shown in 
Figure E.l. The black body temperature can then be determined from the peak 
intensity of an emission line. Since absolute intensity measurements are 
prone to error, fitting the black body curve to several separated peaks of 
optically thick lines is more accurate. The expression for the flux of rad- 
iation I per unit wavelength and unit solid angle is [E.l] 


. 2hc^ l_ 

o “ .5 hc/kTA 

A e 


(E.3) 
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An optically thick line approaches the black body radiation 
curve, I q (T) . 


Figure E.l 



where h Is Planck's constant and c Is the speed of light. The Measurements 
should be performed near the peak of this curve where the sensitivity to 
temperature variation is largest so as to enhance the accuracy. Since at 
high temperature the peak is In the UV and at low temperature it is in the 
IR, it is not always possible to satisfy this condition with one spectrometer. 

Thermal emission allows separate species to be observed since each atom 
or molecule has a characteristic spectrum. The species participating in the 
emission process may be identified by passing the collected light through a 
spectrometer. A measurement of the relative concentrations of these species 
is complicated by the different emission coefficients and optical quenching 
in the optically thin limit and different emissivities and optical depth in 
the optically thick limit. 

The spectral profile of emissions from hot luminous particles may also 
be resolved by using a spectrometer. The temperature of these particles can 
be obtained by fitting the spectrum with the black body distribution given 
in equation (E.3). The particle's temperature may deviate from the gas tem- 
perature particularly in regions of high velocity and high temperature grad- 
ients in the flow. 

A method has been developed to measure hot luminous particle velocities 
using a time of flight technique [E.3]> Two points of known separation in 
the flow are focused onto two different detectors. As the particle passes 
through, its emission is picked up by the detector as two pulses. One can 
determine the velocity by the difference in arrival time between the pulses. 
The spots have to be far enough apart to provide accurate time measurement 
but close enough to contain less than two particles. The vector connecting 
the two sample points in the flow has to be close to the particle trajectory 
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to ensure signal detectability. Problems similar to those associated with 
the classical laser Doppler velocimetry exist for this method, particularly, 
the inability of particles to follow the flow in regions of high accelerations. 

The most serious aspect of thermal emission is its potential as a noise 
source which may obscure other measurements. An estimate of the noise inten- 
sity begins with the radiative transfer equation [E.4] 

>5 


£<Tr + -5r> - \ (,+ -7 V - K * '> 


(E.A) 


where 1^ emission is the intensity, s is the displacement along a ray of 
light, and and are the emission and absorption coefficients per unit 

wavelength respectively. The emission and absorption coefficients may be 
related by invoking detailed balance. For the radiation to be in thermo- 
dynamic equilibrium, the medium has to be optically thick at all frequencies. 
Alternatively, one can assume that collisional processes are much more 
important than radiative. In other words, the excited state has a much 
higher probability of de-excitation by an inelastic collision than by spon- 
taneous emission. A temperature may then be defined that fits the Boltz- 
mann, Saha and Maxwell relations for the excited and ionic populations and 
velocity distributions. The radiative transfer equation at steady state 
may then be rewritten using Kirchhoff's law as [E.4] 


it-SVa’ • 

where I is the black body intensity distribution (E.3) and 

v /i “hc/kTA* 

K A " *A " e > 


(E.5) 


(E.6) 


The equation may be integrated analytically for a homogeneous radiating 
plane layer of finite thickness, d [E.4]. 
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d 

l x (d) - l o 0-e X ) (E.7) 

A schematic of the detection configuration expected on a re-entry vehicle 
is shown in Figure E.2. The hot region between the shock and the wall is 
radiating. It is assumed for simplicity that the gas is homogeneous though 
very sharp gradients exist. This is Justified by the fact that the high 
temperature regions dominate the radiation and will be used for estimating 
the emissions. A lens which is used to collect the signal generated by one 
of the nonintrusive optical techniques also collects background light due to 
thermal emissions. The thermal emissions are a distributed source so light 
is collected from regions other than the focal volume. It may be shown that 
the light is gathered from an equivalent truncated conical volume defined by 
the aperture and the lens as indicated in the figure. Equation (E.7) may 
then be integrated to yield 

l Xa -f l 0 T(F)' 2 (E.8) 

where F is the f number of the collection optics (F- lens focal length/lens 
diameter) and T * K^d is the optical depth (assumed to be «1). The ab- 
sorption coefficient per meter in hot air (T < 20,000°K) is taken from 
Zeldovitz and Raizer [E . 4 ] . 

1 10~ 2 -1 1*0,000 -181 ,000 
K a * 8.6- 10 19 p-T-X 3 e ^ (C Q e T + e T + 

- 158,000 - 169,000 - 108,000 

c 0 e T ♦ C N « T * C N0 e T ) (E.9) 

where p is the density in Kgm/m , T is in °K, X is in meters and C. is the 
concentration of species i. K, and I have been plotted for different 

A O 

temperatures versus wavelength in Figures E.3 and E.l*. The flow conditions 
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Luminous Shuttle 




Figure E.2 Optical configuration showing location of shock, luminous layer. 
Space Shuttle wall with lens and aperture. 
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assumed are taken from the simulation by J.L. Shinn, J.W. Moss and A.L. 
Slmmonds of NASA Langley for a nonaqult Ibrlum flow, fully catalytic wall 
calculation at an altitude of 75 km and at a post Jon along the surface of 
the Shuttle corresponding to 10 nose radii. 

As t example, consider: 
d » 0. 1 m 

X - 510 nm (copper vapor laser) 

F - 10 
T - iO^K 

— Q 

AX » 0.1 x 10 m (state-of-the-art Interference filter) 

From Figures E.3 and E.4, I ■ 2 x iO® W/m^y and ■ 1.1 xIO ^ m The 
optical depth, t. Is 1.1 x 10 1 and equation (E.9) may be used. The 

-3 2 

Intensity radiated Into the aperture Is 2x10 W/m . By comparison, the 

-15 2 

region in the flow at a temperature of 3000° K contributes only 10 W/m 
and may clearly be neglected. 

In summary, thermal emission from gases does not seem to be a promising 
method as far as quantitative measurements are concerned. This conclusion 
arises from the integrated path which makes point observation Impossible and 
the uncertainty regarding the values of the quenching rates and Franck- 
Condon factors which are necessary for unequivocal data reduction. On the 
other hand, thermal emissions contain information which can be utilized 
for nonquant I tat I ve species identification. 

Thermal emissions from gases and particles constitute a noise back- 
ground source for other optical nonintruslve techniques. The first experi- 
ment to be undertaken should include a monitor for thermal emissions. In 
particular, an optical multichannel analyzer would be able to take a complete 



spectrum in a much shorter time than a scanning spectrometer. The spectrum 
would be indicative of the species existing in the flow. Bursts due to hot 
particle; passing through the probe volume can be gated out or directed into 
a different channel for additional processing. 
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APPENDIX F. ELECTRON BEAM FLUORESCENCE 


The electron beam fluorescence technique is used to measure the static, 
vibrational, rotational and translational temperatures of low density gase- 
ous flows as well as determine species concentration, measure velocity and 
perform flow visualization. Reviews of this technique have been written by 
E. P. Muntz [F.l] and this discussion draws heavily on these works. 

A collimated beam of electrons of about 1 mm in diameter is directed 
into a flow, exciting the gas constituents into higher energy levels. The 
species decay after a characteristic lifetime emitting fluorescence. An 
optical detector may be configured to monitor any position along the beam 
to give a spatial resolution of M mm^. The resulting spectrum is charac- 
teristic of the composition, temperature, density and velocity of the gas. 
The many different parameters on which the electron beam fluorescence is 
dependent make it an attractive technique. Its use is usually limited to 
numb' nsities of less than 10^ cm ^ though some investigators have 
extended the range to 10 cm 3 [F.2]. The electron beam loses its colli- 
mation and penetration depth at higher densities and the fluorescence is 
no longer linearly related to density. Ambient temperatures and less have 
been measured and some work at up to 1000°K has also been performed [F.3]. 
The velocity sensitive feature of the electron beam fluorescence technique 
has been applied only in a few cases and then only to hypersonic facilities 
(1000 to 2000 m/sec) [F. 1 *]. 

Nitrogen and air flows have been the most extensively investigated. 

The major emission in both occurs for the first negative system of N^ + at 
391 .1* nm. A 50,000 electron volt beam exciting the first negative system 
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of nitrogen with a I mA/mm of current will produce 10 photons per second. 


Incomplete knowledge with regard to quenching rates and F ran ck- Con don fac- 
tors is the main limiting factor as far as data interpretation is concerned. 

The most straightforward quantitative measurement is the nitrogen 
density. This requires a density independent transfer of energy from the 
electrons to the nitrogen excited state. There are several mechanisms for 
accomplishing this transfer: 1) Collision of nitrogen with primary high 

energy electrons originating from the electron beam; 2) Collision of nitro- 
gen with secondary electrons produced by species ionization due to colli- 
sions with primary electrons; 3) Absorption of electron beam generated 
photons by unexcited nitrogen molecules; k) Energy transfer by collisions 
of nitrogen with other excited particles; 5) Transfer of population to the 
excited state by cascade from higher electronic states. The excited state 
is subsequently depopulated by spontaneous emission and quenching colli- 
sions. The emission is observed as signal; the quenching collisions com- 
pete with emission and reduce the signal. 

If only the primary electrons are assumed to populate the excited state 
then the fluorescence emission, 1^, is related to the gas density, n^, by 
[F.1] 


'lk ' K "g 


Q V 

1+2n 

9 A lk + \ 


(F.l) 


where K is a constant, Q. is the quenching collision cross section. A., is 

99 * x 

the spontaneous emission for the transition of interest, Aj .s the spontan- 
eous emission for all otner transitions and V is the relative velocity 

99 

between the molecules. 

The plot of this relationship is shown in Figure F.l. At low density 
the emission, I ^ ^ , Is essentially proportional to the gas density. At 



F-3 



% 


Figure F.l The Electron Beam Fluorescence intensity saturates at densities 
above 10 16 cm 3 . 
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higher densities the fluorescence intensity is strongly quenched end is 
no longer sensitive to density variations. Host density measurements are 
performed In the linear portion of the response curve. The accuracy of 
density measurement can approach ±1 percent if the calibration procedure 
is carefully executed; lowever, ±5 percent Is more typical. The large num- 
ber of scattered photons reduces photon statistical noise to a low level. 
Inaccuracies are introduced by the quenching process which is both species 
and temoerature dependent. 

Some experimental work has been carried out at higher densities [F.2] 

( 10^ £n lo' ® cm ^) . This is the nonlinear regime and a careful calibration 
is required over the complete range of density and temperature. 

In a nonequilibrium gas, the vibrational aid rotational energy level 
populations may not be ii equi 1 ibrium wi th each other or with the local 
temperature. For example, the vibrational leveis may separately reach an 
equi I ibrium which may be described as a temperature but is different from 
the true translational temperature of the gas. To measure this vibrational 
temperature, fi rst consider a ground state molecule with vibrational level 


v'j' excited to the v 1 vibrational state of the upper electronic level with 
subsequent decay to the v^' vibrational level of the ground electronic state. 
If only spontaneous emission depopulates the upper state and quenching is 
neglected, then the intensity of fluorescence is [F.l] 


I 




q(> ' ,v'J) l n q(vV') 

^ 1 1 V « 1 

vy 1 



(F.2) 


where K is a constant, v is the frequency, n is the number density, anr* q 
is the Franck-Condon factor. The vibrational temperature information is 



which is assumed to follow 


contained explicitly in the number density, n ,,, 

V 1 

a Boltzmann distribution. By monitoring the fluorescence intensity from 
two separable vibrational lines and taking their ratio, one may obtain a 
curve which is temperature sensitive and independent of the constant K. 

The technique has been checked at low densities and up to 300*K with 
a resulting accuracy of ±10% [F.5]. Quenching invalidates equation F.2 at 
higher densities. The number of vibrational levels involved in the calcu- 
lation rises with temperature. This introduces an additional uncertainty 
due to the larger number of Franck Condon factors involved. Higher temper- 
atures also have a bearing on the quenching process. 

Similarly the rotational modes may reach a separate "temperature". The 
expression which provides the rotational temperature is rather involved, 
requi ring an i teration process and will therefore not be given here. Similar 
difficulties to those associated with the vibrational temperature measurement 
are encountered with the rotational temperature measurement and they are 
related to uncertainties in the knowledge of quenching rates and Franck- 
Condon factors. The rotational temperature can in many cases be assumed 
equal to the translational temperature; however this is not always true for 
the vibrational temperature. 

Low density is still a requirement to minimize quenching, but accura- 
cies of better than ±5°K have been achieved at temperatures below 1G00°K [F.6]. 
Large discrepancies exist for temperatures below 300°K. Their origin is 
still unknown. 

The translational temperature is considered to be the true local tem- 
perature. The translational temperature may be obtained by spectrally 
resolving the fluorescence emission of a particular transition. The 
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instrument most frequently used is a Fabry-Perot interfereometer. The line 
shape observed has a Gaussian component which is related to Doppler broad* 
ening due to the random motions of the molecules. The temperature of the 
flow, T, is related to the full width at half maximum, Av^, of the Gaussian 
1 ineshape by [F. 1 ] 


T = 



m 

2k ln2 


(F.3) 


where X is the wavelength of the transition, m is the mass of a single atom 
and k is Boltzmann's constant. Under optimum conditions the technique can 
measure the translational temperature in helium with an accuracy of a few 
degrees K at room temperature. The accuracy in other gases depends on the 
fluorescence intensity and the ratio between the Doppler width and the col* 
lisional and instrumental widths. 

Electron beam fluorescence can be used for velocity measurements by 
observing the Doppler shift of the fluorescence. A comparison of the fluor- 
escence spectrum of a stationary gas with the same transition in moving 
gas gives the Doppler shift, 6Vp. The velocity, U, may then be calculated 
from [F.l»] 

U = 6v D • X (F.A) 


where X is the wavelength of the observed transition. The relative accuracy 
of the measurement depends on the relative size of the Doppler shift, 5Vp, 
in comparison to the combined temperature collisional and instrumental 
broadening, Av. A tracking Fabry-Perot simplifies velocity measurements 
since it eliminates the need to scan the complete spectrum. 

The fluorescence emission also provides a method for density sensitive 
flow visualization. Fluorescence in the close proximity to the electron 
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beam is observed for N^, 0^ and NO. An afterglow is obtained in He flows 
due to the excitation of a metas table state. Hence while the gas fluoresces 
a long way downstream in He flows, in air the electron beam has to pass 
through the position of interest. 

The researchers contacted all agreed that the utility of the electron 
beam technique is limited by insufficient knowledge of quenching rates which 
leads to large uncertainties in temperature measurements. It was felt by 
Professor Larry Talbot of the University of Cal i fomia-Berkeley and Dr. 
William Hunter of NASA Langley that electron beam fluorescence was most use- 
ful as a flow visualization tool. Dr. Robert Cattolica of Sandia-Livermore 
felt that it would be worthwhile to undertake measurements at very high al- 
titudes. At high altitudes (somewhat above those of interest in this study) 
the densities will be low enough so that the data could be interpreted 
quantitatively. At lower altitudes quenching will eliminate this possibil- 
ity, but gross features of the flow may still be observable. 

In summary, the electron beam fluorescence technique is attractive due 
to the many properties which it can measure and the fact that the signal is 
not photon statistics noise limited. For linear operation, it can only be 
used at the highest altitude of interest. At lower altitude one could possi- 
bly apply it to density measurements only after careful calibration. 

Major uncertainties in temperature measurements arise from insufficient 
knowledge of the quenching and Franck-Condon factors particularly at the 
high temperatures existing in the flow. Velocity measurements are best per- 
formed on atoms but have not been demonstrated on nitrogen or other species 
present in the flow. Since a fluorescence signal will also be present, it 
might be useful to use an electron beam at least for flow visual ization 
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and try to extract quantitative data by proper calibration at the higher 
ai titudes. 


References 


F. I E.P. Muntz, Methods of Experimental Physics, 18B . R.J. Emrich editor, 
(Academic Press, New York, 1981). 

E.P. Muntz, AGARDograph 1 32 (1969). 

F.2 J.A. Smith and J.F. Driscoll, AGARD Conference Proceedings No. 193 

on Applications of Non- Intrusive instrumentation in Fluid Flow Research, 
AGARD CP-193 (Harford House, London, 1976), paper 16-1. 

F.3 W.W. Hunter, AIAA J. 8 , 959 0 9 70). 

F.4 V, Natiello, Thesis, Princeton University, 1970 (unpublished). 

M. Becker, F. Robben and R. Cattolica, AIAA J. \ 2 , 1247 0974). 

F.5 G.O. Langstroth, Proc. Roy. Soc. A146 , 166 (1934). 

C.Y. Fan, Phys. Rev. J03, 1740 (195?T. 

J.R. Sheridan and K.C. Clark, Phys. Rev. 140 , 1033 (1965). 

F.6 D. Coe, F. Robben and L. Talbot, in "Rarified Gas Dynamics", R. 

Campargue editor, Abst. 160.C.E.A., Paris, 1978. 

N. V. Karelov, A.K. Rebrov and R.G. Sharafutdinov, in "Rarified Gas 
Dynamics", R. Campargue editor, Abst. 135.C.E.A. , Paris, 1978. 



G-l 


APPENDIX G. RAMAN SCATTERING 


Raman scattering is an inelastic scattering process where the wave 
number shift of the scattered light Is characteristic of the spacing between 
the scattering molecule's rotational, vibrational and/or electronic energy 
levels. Raman scattering can be used to measure the number density of the 
scattering species and the gas temperature. 

For a specific rotational, vibrational and/or electronic energy level 
transition, A + B, the Raman scattered signal strength. In terms of power, is 
given by the following expression 


P R " P o N A R L Q c \ n D 


(G.l) 


where P is the power of the incident laser beam, N. is the number density 

O r\ 

of the scattering species in energy level A, (^) ^ is the differential Raman 
scattering cross-section, L is the length of the incident beam from which 
scattered light is collected, is the solid angle of the collection optics 
H is the efficiency of the collection optics and is the detector quantum 
efficiency. 

The Raman scattered light corresponding to different energy level tran- 
sitions therefore is at different wavelengths and its intensity is propor- 
tional to the population of the initial energy level. This results in a 
spectral shape which is a function of the relative energy level populations 
and which under equilibrium conditions can be related to the gas temperature 
The integrated spectrum in turn is proportional to the total nuihber density 
of the scattering species. 

The sensiti/ity and accuracy of Raman scattering for species density 


measurements was evaluated in terms of the minimum detectable species 
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density for a given signal to noise ratio. The noise was assumed to be 
quantum noise which is equal to the square root of the number of detected 
Raman photons, therefore the signal to noise ratio is given by: 

(S/N) R - (P R At n p /hv) ,/2 (G.2) 

where At is the pulse duration, n p is the number of pulses during the mea- 
surement interval, h is Planck's constant and v is the radiation frequency. 

In order to determine the minimum detectable number density for a given 
signal to noise ratio, the experimental parameters in equations (G.l) and 
(G.2) must be specified. The experimental parameters have been selected to 
give maximum signal strength and to be consistent with any constraints imposed 
by the proposed flow field measurements. The experimental parameters are 
summarized in Table G.l. Two lasers were considered, a KrF laser and a cop- 
per vapor laser. The KrF because of its high average power (l J/pulse at 
10 Hz) and its UV output (249 nm) which results in larger Raman cross-section 
and in greater detector quantum efficiencies. The copper vapor because of 
its high average power (10 mJ/pulse at 5000 Hz). The measurement length L 
is chosen to be 5 mm based on the spatial resolution required to characterize 
the boundary layer profiles which were predicted by the numerical computa- 
tions. It is assumed that for the density measurements the integrated Raman 
spectrum would be used, in which case becomes the total number density of 
the scattering species. Based on the rate of descent of the Space Shuttle, 
it is assumed that an integration time of 10 seconds could be used to make 
the measurements. In Table G.2 are listed the five species which are ex- 
pected in the boundary layer, the Raman transition, the Raman cross-section 
and the minimum detectable density for a signal to noise ratio of 10. Results 
for both the KrF and the copper vapor laser are given. 
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Table G.l 


Values Used for Raman 

Scattering 

Signal to Noise 

Estimates . 



KrF Laser 

Copper Vapor Laser 

Peak Power 

P 

0 

40 MW 

0.33 MW 

Pulse Length 

At 

25 nsec 

30 nsec 

Repetition Rate 

n 

P 

10 2 /sec 

5 x loVsec 

Wavelength 

\ 

249 nm 

510 nm 

Detector Quantum Efficiency 

"o 

0.35 

0.1 

Collector Optics Efficiency 


0.1 

0.1 

Optics Collection Angle 

n c 

0.01 str 

0.01 str 

Scattering Length 

L 

5 mm 

5 mm 



Table G.2 

Minimum Detectable Density by Raman Scattering tompared with Expected Maximum Densities at 75 and 

and 52 km A 1 t i tudes . 

Minimum Detectable Density 





wi th 

S/N- 10 

Maximum 

Density 

Species 

Raman 

Tran si tion 

da /dal * 
dQ/ da|N 2 

KrF 

Copper Vapor 

75 km 

52 km 

N 2 

Vibrational 
2332 cm" 1 

1.0 

1L 

4 x 10 /cc 

1 .6 x 10 ^/cc 

3 x 10 ,6 /cc 

1 x 1 0 1 7 /cc 

°2 

Vibrational 
1557 cm' 1 

1.2 

3.4 x lo’Vcc 

1 . 3 x 10^/cc 

6 x 10 ,5 /cc 

2 x 10 ,6 /cc 

NO 

Vibrational 
1877 cm' 1 

0.55 

111 

7.2 x 10 /cc 

2.8 x 10 16 /cc 

3x 10 ,5 /cc 

1 x 10^/cc 

N 

Not 

Appl i cable 






0 

Electronic 

158 cm 1 , 
226 cm' 1 

1.26 

14 

3.2 x 10 

l.3x 10 16 

6 x 10’5/cc 

2 x 10*^/cc 

* dol 
dfl|N 2 

- 4.4 x I0" 31 cm 2 

/str-molecule 

with 5145^ laser. 





fiCOd dO 



The detectability estimates given in Table G.2 can be compared to the 
calculated species densities at different altitudes, which are also summar- 
ized in Table G.2, in order to determine the altitudes at which Raman scat- 
tering concentration measurements are possible. Based on these results and 
the requirement of a dynamic range of 10.lt can be concluded that with the 
KrF laser N^, 0^, NO and 0 can be measured at all altitudes; whereas with 
the copper vapor laser, can be measured at altitudes up to 60 km while 
0^, NO and 0 can be measured at altitudes up to 52 km. It is important to 
re-emphasize that these detection estimates are made assuming no appreciable 
background noise due to particulates (i.e. scattering or incandescence) or 
thermal emission. 

The calculations predict that the nitrogen mass fraction varies by less 
than 30$ across the boundary layer profile for a given altitude and distance 
along the vehicle, therefore from the point of view of characterizing the 
reacting nature of the boundary layer flow, the nitrogen density is of little 
interest. However, the gas density does change by approximately a factor of 
10 over a given boundary layer profile, therefore a nitrogen density mea- 
surement could be used as a measure of the total gas density. 

For temperature measurements one should use the most populous species, 
if possible, which in this case would be nitrogen. As shown above, using a 
KrF laser the detection of total N^ concentration is possible at all alti- 
tudes. However, if one chooses to resolve the Q-branch spectrum for the 
purpose of curve fitting to determine the temperature, the partitioning of 
the signal among the different vibrational-rotational levels will increase 
the minimum ^ density at which thermometry measurements can be made by a 
factor of 100. In this case thermometry measurements would only be possible 
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at 52 km and below. If a vibrational band ratio or Stokes/anti -Stokes ap- 
proach is used for larger signal strengths, then the partition factor is in 
the worst case 10 and thermometry measurements can be made at all altitudes. 

Assuming temperature measurements could be made from the nitrogen 
vibrational Raman spectrum as discussed above, one would have to be concerned, 
especially near the shock, about non-equi 1 ibrium effects, i.e. whether the 
measured vibrational temperature was in equi 1 i bri urn wi th the gas kinetic 
temperature. 

Another possibility for enhanced signal strengths would be to use reso- 
nance Raman scattering off of 0^ (193 nm) or NO (220 nm) whereby several 
orders of magnitude greater sionel strengths could be realized. As pointed 
out by Dr. Alan Eckbreth of United Technologies, this would require the use 
of a laser with excellent frequency stability. There are, however, no re- 
ports in the literature of reonsnace Raman measurements of 0^ or NO there- 
fore little more can be concluded at this time. 

The concensus of the researchers we interviewed was generally not as 
optimistic as these results, however most ot their concerns were related to 
the potential problem of background noise from particulates and thermal 
emission. There was agreement that characterization of these potential noise 
sources must be done before the final experiment selection and design is 
completed. Because of the weak signals associated with spontaneous Raman 
scattering, there was also agreement that fluorescence and Rayleigh measure- 
ments should be given serious consideration. 
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APPENDIX H. COHERENT ANTI-STOKES RAMAN SCATTERING 


Coherent anti-Stokes Raman scattering (CARS) Is a three photon, in- 
elastic scattering process which is enhanced as the difference in frequency 
between two Incident laser beams nears resonance with an energy level tran- 
sition of the scattering molecule. The scattered radiation is then at the 
Raman anti-Stokes wavelength. The intensity of the CARS signal, 1^, at 
frequency w AS is given by the following equation. 


whe re 


as 


(^) 2 ' 2 * s lxl 2 z 2 n 


f H . 1 ) 


to.. ■ 2u 
AS p 

x - X B + X, 


GO 


S * 


NR 


and 


2Nc ,d Ox 
X R ■ T (“) 


Aw A , 


2 v dfJ / R .5 “ Xo + *Xp 

fi w c Aw-(w - a> c ) Z - i y (oj^ - ) 


Up . IX Wp 

Ip is the intensity of the probe beam at frequency Wp, lj is the intensity 
of the Stokes beam at frequency w c , X I s the third-order nonlinear suscep- 
tibility, Z is the distance over which the beams interact, c is the speed 
of light and n is the combined interaction-collection-detection efficiency. 
The third-order susceptibility is composed of a resonant part Xp and a non- 
resonant part x NR - The resonant part of the third order susceptibility is 
dependent on N, the total '.pecies number density, (^r) R the Raman cross- 
section for the transition characterized by the frequency difference Aw, 
y, the Raman line width, Aj, the fractional population difference between 


the levels in the transition and fi is Planck's constant divided by 2rr . 


The CARS process is nonlinear in incident beam Intensity, number den- 
sity, and interaction length whereas spontaneous Raman scattering is linear 
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In these parameters. Because of this nonlinear behavior, for a given exper- 
imental situation, the CARS signal strength can be larger or smaller than 
the corresponding spontaneous Raman signal depending on the number density. 
The CARS process must also satisfy a phase matching requirement (l.e. con- 
servation of momentum). Phase matching is possible with one-sided optical 
access, however a wide angle crossed beam configuration must be used. This 
results in a shortened interaction length and limited traversing capability. 

In order to estimate the CARS signal strength, one can assume that the 
Stokes and anti-Stokes beams have the same diameter. Equation (H.l) can 
then be written in terms of the CARS power 

'as • l*l 2 'p " s * 2 1 <"- 2 > 

C 

There are two limiting cases to consider when estimating the minimum 
detection limits. One case is the detection of a major species for which 
X d >Xud* and the other is the detection of a minor species for which 
Xd < XkiD * We will consider the measurements to be quantum noise limited in 
both cases. In the former case the greatest source of noise is the quantum 
noise in the signal itself, therefore 

P AC At n ... 

< S '">CARS ' (- 5 V- E) <H - 3) 

where J X j 2 “ Xj^ + X'^ 

In the latter case the major noise source is the quantum noise in the non- 
resonant background. Therefore the signal is the contribution due to 

2 

Z *R*NR an< * t>1e no * se * s th * squ ar « root of the contribution due to x NR * 

2 2 

The ratio of these two reduces to equation (H . 3) where |xP * **Xr * s used 
in equation (H.2). Therefore the difference in the detection limit for the 
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major species and the minor species cases differs only by approximately a 
factor of two. 

The minimum detectable nitrogen density for a signat to noise ratio of 

10 can now be determined using equations (H.2) and (H.3)* The assumed 

experimental parameters are listed in Table H.l. The resultant detectabil- 
1 5 3 

ity limit is '.D /cm . This esti..iate is for 1500°K, however away from the 
surface the .emperatures are expected to be much higher. As the temperature 
increases, the rotational population spreads out over many more levels which 
causes a reduction in the CARS strength through the Aj term in the third 
order susceptibility. And as Eckbreth pointed out, the contribution of ad- 
jacent lines will also be less at high temperature which will further reduce 
the CARS strength. 

If we compare this result to the calculated densities of N^, 0^, NO and 
0 at different altitudes we conclude that N^ can be det* ted at all altitudes 
and that 0^, NO and 0 can be detected at 52 kir. and below. Again, the sensi- 
tivity will be less in the hotter reaions away from the surface. In order 
to relate the CARS signal strength to the density, ! t it necessary to know 
the gas temperature, unless a particular transition can be selected for which 
the population difference term is temperature independent. Determination of 

the temperature from the N 2 CARS spectrum requires resolving the spectrum 

1 7 3 

which increases the minimum detection limits to 10 /cm because of reduced 
away from the peak and the necessity to use broad band CARS. Thermometry 
measurements therefore will only be possible ;t the lowest altitude. 

Because the detection limit analysis predicts adequate sensitivity at 
only the lowest altitudes and because of the severe limitation due to one 
side optical access, it is concluded that CARS is not viable for the proposed 



re-entry vehicle measurement. Several of the researchers interviewed have 
had extensive experience using CARS diagnostics and their comments were in 
general agreement with this conclusion. 


Table H.i 

Values Used for CARS Si gnat to Noise Estimates. 


First laser wavelength 
Second laser wavelength 
Observed wavelength 
First laser peak power 
Second laser peak power 
Nonlinear susceptibility 
Combined efficiency 
Beam diameter 
Pulse duration 

Total number of pulses integrated 
Interaction length 


X - 

P 

X « 

s 



P - 

P 

P » 

s 

Xp “ 3*2x10 cmVe r g at 1 atm, 

n * 

D - 

At - 

n * 
P 

2 “ 


532 nm 
607 nm 
480 nm 
10 MW 
1 MW 
1500°K 
0.01 
200 ym 
10 nsec 
100 

200 ym 
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APPENDIX I. LASER DOPPLER VELOCIMETER 


The laser Doppler velocimeter Is useful for measui Ing fluid velocities. 
Laser light scattered by a particle moving with the flow experiences a Dop- 
pler frequency shift proportional to the velocity of the particle. A mea- 
surement of the shift in frequency then gives the particle velocity. The 
most commonly employed configuration uses two laser beams crossed in the 
sample region to produce interference fringes (Figure l.l). As the particle 
traverses this region it crosses the fringes (Figure 1.2) at a rate deter- 
mined by the veloci ty, resul ting in scattered light which is modulated at 
the crossing frequency. Spatial resolution Is defined by either the beam 
crossing volume or the collection optics with submillimeter resolution easi- 
ly obtainable. A second approach employing a Fabry-Perot interferometer may 
be useful at high fluid velocities. The laser Doppler velocimeter is non- 
intrusive and insensitive to other fluctuating flow parameters such as 
density and temperature. Velocity ranges from 10 ^ cm/sec to 10^ m/sec have 
been measured. The difficulties associated with the use of the laser Doppler 
velocimeter are primarily due to the fact that the flow has to be seeded and 
that the particles do not track, the flow in regions of high acceleration or 
low density. 


The detected signal arising from light scattered by the particles as 
they move through the set of interference fringes may be found from equation 


1.1 [ 1 . 1 ] 


\ 

o 

2 sin a/2 


(l.l) 


Here Is the particle velocity In the plane of the two laser beams and 
perpendicular to their bisector, o is the intersection angle, X Is the 
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laser wavelength, and Vg is the beat frequency produced. The distance be- 
tween the fringes is x^. The product of the output of the frequency coun- 
ter, Vg, and the fringe spacing, x^, equals the velocity. It is clear that 
for a given laser frequency and beam intersecting angle, the beat frequency 
is directly proportional to velocity. Commercial laser Doppler velocimeter 
units can handle counting rates up to 200 MHz. State-of-the-art counters 
can go as high a 1 GHz which is also approximately the largest bandwidth of 
available photomultipliers. Counter availability is therefore the limiting 
factor as far as high velocity measurements are concerned. Table 1.1 lists 
the counting frequencies, Vg, for different beam intersecting angles, a, 

for an argon laser operating on the A = 51^.5 nm line and a particle velo- 

o 

city of U. =5000 m/sec. 

Table 1.1 

Fringe Spacing and Counting Frequencies for .51*»5y Laser Beams 
Intersecting at Angle at ana 5000 m/sec Particles. 


q° 

Xf (ym) 

VgjHH 

*♦5 

0.7 

7000 

20 

1.5 

3333 

5 

6 

835 

1 

30 

166 


From the table it is clear that the smaller the angle the more manageable 
the counting rate. If this method is to be used with a re-entry vehicle, 
very small angles will be required. An alternative approach to reduce the 
counting rate is to shift the frequency of one of the laser beams so the 
fringes scan in the direction of motion. This eliminates part of the Doppler 
shift and is useful only if a narrow range of velocities are to be measured. 
A large value of interaction length permits one to focus several detectors 
at different positions along the interacting volume so as to obtain a velo- 
city profile. Clearly the shallower the angle o the longer the interaction 
length becomes. 
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The second approach circumvents the high counting rate problem by re- 
solving the scattered light spectrally. A Fabry-Perot etalon is inserted 
between the intersecting volume and the detector. Two peaks are registered 
by the detector as the etalon is scanned as shown in Figure 1.3. They cor- 
respond to the Doppler downshifted and upshifted light scattered by the two 
laser beams off the moving particles. The frequency difference between the 
two peaks is v„ in equation (1.1). Av is the free spectral range of the 
Fabry-Perot which provides a frequency caliDration. While an upper bound 
limits the velocity measureable with the optical mixing conf igurat ion, a 
lower bound exists if the scattered light is to be spectrally resolved. 

When both scattered signals are put through the same etalon the minimum mea- 
sured velocity is limited by the Rayleigh criterion: a shift is measureable 

only if it is larger than the width of the spectral line, Av^v^, otherwise 
the two peaks coalesce. To overcome the limitation, a reference beam config- 
uration may be employed whereby just one beam interrogates the flow while 
the other beam is used as a frequency reference. The scattered light from 
the flow and the reference beam are now directed into different etalons so 
that the peaks cannot overlap at small velocities. A single etalon may be 
employed if the frequency of the reference beam is shifted using a Bragg cell. 

Two additional techniques which have been used and are related to laser 
Doppler velocimetry should be mentioned. One is the time of flight veloci- 
meter [1.2] which may be used if the flow direction is known. Two parallel 
laser beams are directed into the flow. A detector monitors the light pulses 
scattered off the particles as they pass through the beams. The velocity is 
deduced from the known separation of the beams and the time difference be- 
tween the light pulses. 
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In a related technique discussed by Or. Don Holve of Sandia Labora- 
tories the pulse shape of a particle traversing a single laser beam Is 
resolved temporally and is independent of the particle size or trajectory. 
The time between two fractional values of the peak amplitude is related to 
the particle's velocity. The intensity of the pulses can also be used to 
obtain the particle size distribution, however velocity-size correlations 
cannot be determined [ 1 . 31 - 

Two extensions of the LDV technique were mentioned in our meeting at 
Sandia. In the first a holographic lens is employed to prolong the fringe 
pattern for multiple joint measurements. The second is a time of 
flight measurement which takes advantage of gas breakdown by laser focusing. 
The emitting plasma follows the flow and is monitored by two detectors which 
furnish the elapsed time. 

The seeding problem poses the greatest difficulty as far as applying 
the LDV to measurements on board a re-entry vehicle. Particles injected 
through the wall of the craft would tend to stay close to the wall. Natu- 
rally occurring seed however may be present or emitted from the ceramic 
cover of the re-entry vehicle and these may be used for LDV measurements. 

The interference fringe detection schemes are more prevalent and easier to 
use. A very narrow angle a will be required to reduce the count frequency, 
but this also has the advantage of permitting multiple point measurements. 

Thermally excited fluorescence emission and particulate luminosity may 
mask the signal. For these reasons it seems advisable as a first step to 
monitor fluorescence and particulate emission. This may be augmented by a 
laser to determine whether particles exist in the flow. Only after the en- 
vironmental effects have been documented should an LDV measurement be at- 


temped. 
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Figure 1.1 A Laser Doppler Velocimeter heterodyne detection configuration. 



Figure 1.3 Fabry Perot output looking at light scattered off particles 
illuminated by two laser beams. 
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APPFNDIX J. RESONANT DOPPLER VELOCIMETER 

The Resonant Doppler Velocimeter (RDV) [J . 1 ] has been used for visual- 
ization and for quantitative measurements of velocity, temperature and den- 
sity of a gaseous flow. These are achieved by shining a single frequency 
laser beam into a flow which is seeded with an atomic or molecular species. 
The laser is tuned through the absorption frequency of the seeded species. 

An optical detector may be configured to monitor fluorescence from any posi- 
tion along the beam to provide good spatial resolution. The velocity is 
determined by observing the Doppler shift of the absorption frequency. Spec- 
troscopic absorption line broadening mechanisms furnish information regarding 
the static temperature and density of the moving gas. 

The technique has been demonstrated in supersonic and hypersonic flows. 
Unlike particles needed for the standard laser Doppler velocimeter, tracer 
atoms injected into the flow follow it closely. But the flow still has to 
be seeded, a fact which could raise experimental difficulties. 

The velocity of the flow is measured by the Doppler shift Av^ of the 
absorption line. 

Av D =V o*c (J ‘ 1) 

A 

where V q is the absorption frequency in the atomic rest frame, k is a unit 
vector in the direction of the laser beam, U is the flow velocity vector and 
c is the speed of light. The method measures the velocity vector in the 
direction of the laser beam. A reference has to be furnished to be able to 
measure the absorption frequency in the atomic rest frame. This can be pro- 
vided by performing spectroscopy in a stationary gas in a cell or in an 


atomic beam device. 
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The temperature and density of the flow give rise to frequency broad- 
ening of the spectral line. The temperature effect results In a lineshape 
which is a Gaussian and the collision frequency, which is directly propor- 
tional to dens i ty, causes a Lorentzian lineshape. Their combination is a 
convolution which is known as the Voigt profile. The observed profile may 
be deconvolved to furnish the separate broadening components. The tempera 
ture, T, of the flow is given by 


t / AV G V M 
T " ( 2 * 2k ln2 


(J.2) 


where Av_ is the full width at half maximum of the Gau-s-an lineshape, X 
G 

and M are the absorption wavelength and mass of the particular species and 
k is Boltzmann's constant. The density of the gas, n, is given by 

IT Av, 


n = 




o u 

where Av^ is the full width at half maximum of the Lorentzian lineshape, 0 
is the collision cross section and u is the relative velocity in the center 
of mass system of the tracer atom and its colliding partner, a is tempera- 
ture independent if the hard sphere potential is used. In practice one has 
to find the temperature dependence of o by performing a calibration [J.2]. 

The resonant Doppler velocimeter may also be used for flow field visual- 
ization. Regions of higher density in the fiow fluoresce more bright. In 
addition, by tuning the laser to a particular Doppler shifted absorption 
frequency, one can highlight a particular velocity field. 

The many properties which can be measured by using the resonant Doppler 
velocimeter make it an .attractive technique. The major difficulty lies in 
identifying the proper tracer atom and seeding it into the flow. Sodium 
atoms have been used in nitrogen and helium flows but sodium would tend to 
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react with the oxygen existing in the flow field surrounding the vehicle. 
Sodium does occur naturally in the atmosphere. At the highest altitude of 
interest (80 km) its concentration ov 10 cm is four orders of r..-agnltude 
below that used in the laboratory. At lower altitudes its concentration is 
negligible and the fluorescence signal will be masked by Rayleigh scattering. 

Laser induced fluorescence has been demonstrated on N, 0 and NO sug- 

gesting that they could be used as tracers for applying the resonant Doppler 
velocimeter. The major difficulty is associated with producing a tunable 
laser system with a narrow linewidth and an output frequency which matches 
one of the ultraviolet absorption frequencies of the gases present. 

References 

J.l M. Zimmermann and R.B. Miles, Appl . Phys. Lett. 32, 885 (1981). 

J.2 M. Zimmermann and R.B. Miles, J. Phys. B: At. Mol. Phys. 14, L85 (1981). 



(u',(/,f)u'i(*.i) 

* autocorrelation in the near wall region defined 

/ 

by Eq. 19 

(uic'Xii) 

m turbulent flux of a passive additive 

«c(*l) 

“ characteristic velocity appearing ; the tur- 
bulent model for (eric'); u c * L(- / T(xi ) 

X 

m coordinate vector with com pone, .s sq, xg, xj 

+ 

m superscript denoting a quantity made dimen 
sionless using the wall parameters u* and v 

Qrook LaMar* 

a 

m defined by Eq. 23 

«c 

“ characteristic size of the concentration layer 
near the mass transfer interface (Eq. 7) 

V 

“ kinematic viscosity 

th 

■ characteristic relaxation time for the velocity 
autocorrelation in a frame of reference moving 
with the average velocity 

(tm) 

* mean period between burst 

a 

“ semi-infinite spatial domain 


LITERATURE CITED 

Abromowitz. M . and I A. Stegun, Handbook of Mathematical Fund tom, 
NBS, App. Math Serin, 55 (1964) 

Berman, N. S , "The Effect of Sample Probe Size on Sublayer Period in 
Turbulent Boundary 1 -ayers," Chem. Eng. Comm., 5,337 (1980). 

Blackwelder, R. F., “The Bursting Process in Turbulent Boundary Layers," 
Coherent Structure of Turbulent Boundary Layers, eds., C. R. Smith 
and D. E. Abbott, AFOSR, Lehigh Univ.. 211 (1978). 

Campbell, J A. and T. J. Hanratty, "A New Theory for Turbulent Mass 
Transfer at a Solid Boundary," submi.ted to A. I Ch.E. }. (1981) 

Cantwell, B J., "Organized Motion in Turbulent Flow," Ann. Rev. Fluid 
Mech., 13, 157 (1981). 

Coles, D., “A Model fur Flow in the Viscous Sublayer," Coherent Structure 


of Turbulent Boundary Layers, eds., C R. Smith and D. E. Abbott, 
AFOSR, Lehigh Univ , 462 (1978). 

Kreplin. H. P., and H. Eckelmann. "Behavior of the Three Fluctuating 
Velocity Components in the Wall Region of a Turbulent Channel Flow,” 
Phys. Fluids, 22 (7). 1233 (1979) 

Kreplin. H. P., and H Eckelmann, "Propagation of Perturbations in the 
Viscous Sublayer and Adjacent Wall Region," J. Fluid Mech. , 95, 305 
(1979b). 

MacLeod, N., and J W Por.ton, “A Model lor Turbulent Transfer Processes 
at a Solid-Fluid Boundary, " Chem. Eng Set , 32, 483 (1977) 

Monin, A. S., and A. M Yaglom, Statistical Fluid Mechanics: Mechanics 
of Turbulence, 1, MIT Press, Cambridge (1971) 

Nakagawa, H., and I. Nezu, “Bursting Phenomenon Near the Wall in 
Open-Channel Flows and its Simple Mathematical Model,” Memoirs, 
Foe. Eng.. Kyoto Univ.. 40, 213 (1978). 

Petty, C. A., “A Statistical Theory for Mass Transfer Near Interfaces," 
Chem. Eng Set., 30, 413 (1975) 

Petty, C. A., and P. E. Wood, “The Effect of Turbulent Mixing on Mass 
Transfer Near Fluid-Fluid Interfaces for Large Schmidt Numbers," 
Chem. Eng. Comm , 7, 261 (1980a) 

Petty, C. A., and P E. Wood, "The Effect of Turbulent Mixing on Physi- 
cochemical Absorption Near Interfaces," Chem. Eng. Comm., 7, 275 
(1980b). 

Shaw, D. S., and T. J Hanratty, "Turbulent Mass Transfer Rates to a Wall 
for Large Schmidt Numbers," AfChE J., 23, 28 (1977a) 

Shaw, D. S , and T. J. Hanratty, "Influence of Schmidt Number on the 
Fluctuations of Turbulent Mass Transfer to a Wall,” AlChE J , 23, 160 
(1977b). 

Sideman, S„ and W V. Pinczewski. "Turbulent Heat and Mass Transfer 
at Interfaces: Transport Models and Mechanisms," Toptcs in Transport 
Phenomena, C. Cutfinger, ed , Hemisphere, 47 (1975). 

van Dongen, F. C., A. C. M Belhaars, «nd D A DeVries, "A Periodic In- 
termittent Model for the Wall Region of a Turbulent Boundary l^yer." 
fnl. ]. Heat Mass Transfer, 21, 1099 (1978) 

Yao, H T.CP. Chen. R W. Snellenberger, P E Wood, and C A Petty, 
"The Effect of a Unite Propagation Velocity on Turbulent Mass Transfer 
Near a Rigid Interface," Chem. Eng. Comm.. 12 (1981). 

Manuscript received July 27, 1981, revision received November 2. and accepted 

Nov em b er 5 . 1981 




Evaluation of a Stochastic Model of Particle Dispersion in a Turbulent 
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The accumulation and correlation of data for empirical com- 
putations of turbulent particle dispersion is complicated since 
dispersion is influenced by both particle and turbulence properties. 
This difficulty can be circumvented by computing dispersion di- 
rectly, using a stochastic particle dispersion model The stochastic 
method requires an estimate of the mean and turbulent properties 
of the continuous phase. Particle trajectories are then computed 
using random sampling to determine the instantaneous properties 
of the continuous phase, similar to a random walk calculation Mean 
dispersion properties are obtained by averaging over a statistically 
significant number of particle trajectories. 

Several stochastic particle dispersion models have been proposed 
Yuu et al. (1978) use a stochastic dispersion model, which employs 
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empirical correlations of mean and turbulent properties, to analyze 
their measurements of particle dispersion in jets. Gosman and 
loannides (1981) describe a more comprehensive approach, pre- 
dicting both flow properties (using a k-t model) and dispersion The 
latter procedure is attractive since k-( models vi.'ld satisfactory 
predictions for many of the boundary layer type flows that are 
encountered with dispersion problems 

The objective of the present investigation was to reexamine the 
data of Yuu et al. (1978) for particle dispersion in air jets using a 
stochastic dispersion model similar to Gosman and Ioanmides 
(1981). Shearer et al. (1979) have demonstrated that the present 
k-t model provides good predictions of existing measurements of 
mean and turbulent properties within jets in the following, the 
present stochastic dispersion model is also calibrated by comparison 
to predictions by Hinze (1975) and measurements by Snyder and 
Lumley (1971 ), for simpler flows, prior to applying the model to 
the jet dispersion data of Yuu et al. (1978). 
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Continuous Phase 

Computations of the properties of the continuous phase were 
only required for the tests of Yuu et al. (1978). The jet dispersion 
experiments involve a particle laden round air jet in;- cted into still 
air which can be modeled as a steady axisymmetric boundary layer 
flow. Particle mass loadings in the jet were 0. 1 -0.4% which is suf- 
ficiently small so that the particles had a negligible effect on mean 
and turbulent gas-phase properties. Since nozzle exist Mach 
numbers were less than 0.3, density variations, kinetic energy and 
viscous dissipation were neglected with little error. Jet Reynolds 
numbers were 9,500-48,000 which implies that molecular transport 
can be ignored in comparison to turbulent transport 

Under these assumptions, the governing equations for the con- 
tinuous phure, including model transport equations for k and t, 
are: 

~ + ™M-D(u)-D(C)-0 (1) 

dx r dr 

D(k) - Mr - pt (2) 

D<<l ' c '“"(i)(^) ! - c ' ! T (3) 


Particle Dispersion 

Following the method of Cosman and loannides (1981) the 
motion of the particles is tracked as they interact with a succession 
of turbulent eddies, each of which is assumed to have constant flow 
properties. Velocity fluctuations were assumed to be isotropic with 
a Gaussian probability density distribution having a standard de- 
viation of (2k/3) ,/2 The local distribution is randomly sampled 
when a particle enters an eddy to obtain the instantaneous velocity 
asu* = u + u'. 

A particle is assumed to interact with an eddy for a time which 
is the minimum of either the eddy lifetime or the transit time re- 
quired for the particle to cron the eddy. These times are estimated 
by assuming that the characteristic size of an eddy is the dissipation 
length scale, similar to Gosman and loannides (1981) 

L, » Cl lA k 3 ' 2 /( (7) 

Gosman and loannides ( 1981 ) compute the eddy lifetime as t, * 
L,/|u'|. however, we found better agreement with measurements 
by departing from this particle and employing 

If = L,/(2k/3) l/2 (8) 

The transit time of a particle was found using the linearized 
equation of motion fo, a particle in a uniform flow 

h * — r In (1 - L,/(t|u* — u p |)) (9) 


where 


DM 


d<t> d4> 

pu~ + pv — 
dx dr 


_1 _d 
r dr 



(4) 


and all dependent variables denote time-averaged quantities. The 
specification of the boundary conditions for Eqs. 1-4 and their 
numerical solution is described by Shearer et al ( 1979). Solution 
of these equations yields u,v.C,k and ( throughout the flow 
field 

A conservation equation for particle concentration appears in 
Eq. 1 , assui, ng that the particles have the same local mean velocity 
and turbulent diffusivitv as the gas phase. This corresponds to the 
locally homogeneous flow (LHF) approximation for two-phase 
flows considered by Shearer et al. (1979). LHF calculations were 
considered since this approximation is frequently invoked. Fur- 
thermore, comparing actual dispersion with the LHF limit provides 
an indication of the effect of the inertial properties of the particles 
on particle dispersion. 


where ti* — u p is the velocity at the start of the interaction. When 
L,> t\u” - u p |. the linearized stopping distance of the particle 
is smaller than the characteristic length scale of the eddy and Eq. 
9 has no solution In this case, the eddy has captured the particle 
and the interaction time is the eddy lifetime. 


CALIBRATION OF MODEL 

Since the stochastic mode! involves some rather arbitrary se- 
lections of length and time scales, it was calibrated similar to 
Gosman and loannides (1981) using the fundamental dispersion 
results of Hinze (1975) and Snyder and Lutnley (1971). Assuming 
constant turbulent diffusivity, Hinze (1975) developed an analytical 
expression for the diffusion of "marked" fluid particles introduced 
at a constant rate from a point source into a homogeneous isotropic 
flaw. Comparable results for the stochastic model were obtained 
by fixing u,k and <. The corresponding turbulent diffusivity of the 
analytical expression was taken as C„pk 2 /t, which is consistent with 


Particl* Motion 

Particle trajectories were determined using a Lagrangian for- 
mulation of the governing equations. Since p,,/p > 200, for the 
measurements to be considered, it is reasonable to neglect virtual 
mass, Bassett forces, Magnus forces, etc Under these assumptions, 
conservation of momentum yields 

+ 151 

where particle motion is based on tlw instantaneous velocity of the 
continuous phase. Gravitational forces were negligible for condi- 
tions examined by Hinze (1975) and Yuu et al (1979), but are in- 
cluded since the effect was appreciable for the experiments of 
Snyder and Lumley (1971), The position cf the particle was de- 
termined from 

m <V * = 1.3 («) 

Equations 5 and 6 were numerically integrated using a second- 
order Runge-Kutta algorithm The trajectories of at least 1,000 
particles were computed and averaged to obtain dispersion prop- 
erties. 



Figure 1. Analytical and stochastic solutions lor dispersion of smalt particles 
In homogeneous Isotropic flow with long diffusion times. 
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Figure 2. Predicted and me Mured particle diaper tlon in a uniform, grtd-gun- 
aratad turbulent now. 


the original derivation. The particles were assumed to be small, to 
correspond with fluid particles. Therefore, the interaction time is 
always I, and particle motion is identical to the fluid motion in each 
eddy. 

The radial variation of particle concentration downstream of 
the point source is illustrated in Figure 1 for the conditions exam- 
ined by Gosman and Ioannides (1981). Predictions for both sto- 
chastic dispersion models are shown, along with the analytical result 
(C c was taken from the analytical expression in all cases in order 
to obtain an absolute comparison). The present stochastic predic- 
tions agree favorably with the analytical expression, while the 
original proposal of Gosman and Ioannides (1981) tends to un- 
derestimate the rate of dispersion. (It should be noted that it was 
not possible to reproduce Gosman and Ioannides ( 1981 ) computed 
results for this flow.) 

The second calibration employed the measurements of Snyder 
and Lumley (1971). These experiments involved the dispersion of 
individual particles which were isokinetically injected into a uni- 
form turbulent flow downstream of a grid. In this case, interaction 
times involved both t, and t,. The comparison between predictions 
and measurements is illustrated in Figure 2 The stochastic models 
provide encouraging agreement with the measurements. The two 
models yield nearly the same results for heavier particles, where 
interaction times are mostly fixed by t,, since this aspect of both 
models is the same. The present model yields somewhat improved 
predictions for light particles, where t, largely controls the inter- 
action time. 


PARTICLE LADEN ROUND JET 

The experimerts of Yuu et al. (1978) involved an air jet con- 
taining nearly monodisperse fly ash particles injected into still air. 
The nozzle was shaped according to the specifications of Smith and 
Wang (1944) to yield a uniform outlet velocity. Although Yuu et 
al. (1978) assumed that gas and particle velocities were identical 
at the nozzle exit, particle trajectory calculations for the stated 
nozzle shape indicated that this was not the case. Therefore. |iarticle 
velocities ?! the nozzle exit were computed to define initial con- 
ditions for the present stochastic dispersion calculations This in- 
volved integrating Eq 5 and 6 assuming that particle and gas ve- 
locities were identical at the nozzle inlet Gas velocities were de- 
termined as a function of distance through the nozzle, assuming 
frictionless flow for the given nozzle shape ?nd neglecting the small 
contribution of particle drag on the gas (low properties Compu- 
tations downstream of the ikjzzIc exit then proceeded as described 
earlier. 

Representative predicted and measured profiles of mean axial 
gas velocities and particle concentrations are illustrated in Figure 
3 Tlie measurements were grouped into ranges of i id. therefore, 
predictions are shown for the limits of these ranges Data scatter 
has been smoothes! by only showing mean values for the mea- 
surements. 

The k-< flow msxlel > lelds fair predictions of mean gas velocities, 
particularly in view of the scatter of tlie data Better perfnrmaner 
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Figure 3. PresHctaC and measured axial gas velocitlM and particle cone en- 
list Ions in a dust laden air jet with particle diameter of 20 pm and particle 
density of 2,000 kg/m*. 


of this flow model was generally observed by Shearer et al. (1979) 
for constant density jets. The stochastic dispersion model, using the 
computed initial particle velocity, yields quite good agreement 
with the particle concentration measurements. Similar computa- 
tions assuming equal particle and gas velocities at the nozzle exit 
resulted in narrower concentration profiles indicating a significant 
influence of initial particle conditions on measured dispersion 
properties. The predictions of the LHF model overestimate the 
dispersion of particles, therefore, particle inertia is important for 
these test conditions The particle trajectory computations also 
yielded significant differences in particle and gas velocities re- 
sulting in maximum particle Reynolds numbers on the order of 
50 


CONCLUSIONS 

A stochastic model of particle disjxrsion by turbulence, proposed 
by Gosman and Ioannides (1981), has been evaluated The method 
employs a k- r model to estimate turbulence properties Disjiersion 
is determined by computing particle motion, with random sam- 
pling to obtain instantaneous flow properties, for a statistically 
significant number of particle trajectories 

The stochastic model yielded good .esulls, particularly when 
eddy lifetimes were evaluated using Eq 8 The advantages of the 
stochastic method are that effects of large relative velocities lie- 
tween the particles and the flow, drag properties at Reynolds 
munliers greater than the Stokes flow regime, and the variation of 
local turbulence properties can be readily handled — at least for 
IxMindary layer flows. 

Additional evaluation of the stochastic model would lx- desirable 
torrent [^rescript ions for eddy length and lime scales are somew hat 
arbitrary and further refinement is needed The use of a linear drag 
law to estimate I, is suspect for the high [article Reynolds munliers 
often encountered in |iractiee Finally, while the assumption of 
isotropic turbulence in the stochastic model is consistent with the 
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k-t turbulence model, this approximation ir questionable Jot most 
turbulent flows of interest (or dtspetMor, problems 
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NOTATION 

C = mean particle concentration 

Ci , = drag coefficient -24U+ Be 1 / 6), He. fir < 1 .000; = 

0.44. fir > 1 C|fi 

C, ~ turbulence mode? constants C u = 0.OT, C,| = 1.44. 
C.z = , 39 

d ~ nozzle diameter 

d p = puftid" diameter 

g = acceleration of gravity 

k = turbulence kinetic energy 

L r = Jissijiatton length scale 

r = radial distance 

fir — Reyp f, t !$ number = p\u" - 'i p \d v p 

t = time 

t, - eddy lifetime 

t, = transit time of particle through eddy 

u = mean axial gas velocity 

u, .Upt = gas and particle velocity c npuients 

r- = mean radial gas velocity 

% .ti = axai distance, component of distance 

y = vertical distance 

Greek Letters 

t = dissipation rate of turbulence kinetic criergv 

a = visCAsitv 


p, = turbulent vises * = Cflfiit 2 ? * 

p = density 

- turbuiciit Prandtl Schmidt numbers: cr„, = <U = I, a, 
= 13. <r< =07 

r = particle relaxation time = 4p ? <f P / : (3fi£-n|u - u p |) 

0 ~ generic property 

Subscripts 

i, = centerline value 

t = ompon»nt in direction i 

p = particle property 

o - initial condition 

Superscripts 

(“) = vector quantity 

(-1 = mean quantity 

(') = fluctuation with respect to mean quantity 

{') = instantaneous value 
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